Dynamics in dense hard-sphere colloidal suspensions 
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The dynamic behavior of a hard-sphere colloidal suspension was studied by X-ray Photon Correla- 
tion Spectroscopy and Small Angle X-ray Scattering over a wide range of particle volume fractions. 
The short-time mobility of the particles was found to be smaller than that of free particles even 
at relatively low concentrations, showing the importance of indirect hydrodynamic interactions. 
Hydrodynamic functions were derived from the data and for moderate particle volume fractions 
($ < 0.40) there is a good agreement with earlier many-body theory calculations by Beenakker 
and Mazur [C.W.J. Beenakker and P. Mazur, Physica A 120, 349 (1984)]. Important discrepancies 
appear at higher concentrations, above $ « 0.40, where the hydrodynamic effects are overestimated 
by the Beenakker-Mazur theory, but predicted accurately by an accelerated Stokesian dynamics 
algorithm developed by Banchio and Brady [A.J. Banchio and J. F. Brady, J. Chem. Phys. 118, 
10323 (2003)]. For the relaxation rates, good agreement was also found between the experimental 
data and a scaling form predicted by Mode Coupling Theory. In the high concentration range, 
with the fluid suspensions approaching the glass transition, the long-time diffusion coefficient was 
compared with the short-time collective diffusion coefficient to verify a scaling relation previously 
proposed by Segre and Pusey [P.N. Segre and P.N. Pusey, Phys. Rev. Lett. 77, 771 (1996)]. We 
discuss our results in view of previous experimental attempts to validate this scaling law [L. Lurio 
et al, Phys. Rev. Lett. 84, 785 (2000)]. 

PACS numbers: 83.80.Hj, 61.05.cf, 64.70.pv, 64.70.qj 



I. INTRODUCTION 

The dynamical behavior of colloidal suspensions is a 
very rich research area with borderlines to many fields 
of fundamental research and important industrial appli- 
cations. High density colloidal suspensions provide, for 
instance, invaluable model systems for the study and un- 
derstanding of dynamics in atomic glasses. Prototyp- 
ical model systems consist in suspensions of spherical 
particles with low size polydispersity. They are stabi- 
lized against aggregation due to van der Waals attractive 
forces by coating the surface with a short-chained poly- 
mer (steric stabilization) or with a charged ionic layer 
(charge stabilization). The present study focusses on the 
dynamics of sterically-stabilizcd suspensions. For such 
suspensions, the inter-particle forces are well described 
by a hard-sphere interaction potential with no detectable 
long-range interactions and an infinite repulsion when 
two particle centers are separated by one diameter. 

The experimental study of dynamics in dense colloidal 
suspensions was pioneered byP. N. Pusey, W. van Megen 
and collaborators (see e.g. [l|-@), using Dynamic Light 
Scattering (DLS). The phase behavior of a hard-sphere 
suspension depends on a single parameter - the packing 
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fraction, or particle volume fraction, In the low vol- 
ume fraction limit, with $ on the order of a few percent 
or less, the dynamics of individual particles is essentially 
Brownian. The relaxation times measured by DLS yield 
a (J- independent diffusion coefficient Dq equal to that of 
free particles i.e. the Stokes-Einstein free diffusion co- 
efficient. However, as soon as the volume fraction is 
increased, the dynamics is slowed down by both direct 
interactions between the particles and by indirect hydro- 
dynamic interactions mediated by the solvent. These in- 
teractions are highly dependent on the structural proper- 
ties of the system and thus on the scattering vector q. In 
order to minimize effects introduced by multiple scatter- 
ing of light, the DLS studies use elaborate refraction in- 
dex matching procedures between the colloidal particles 
and the solvent or complex scattering techniques such as 
two-color DLS (TCDLS) 0. The TCDLS work described 
in Ref. Q studies the dynamics of hard-sphere colloidal 
suspensions. The q- and time-dependent short-time diffu- 
sion coefficient Ds{q,t) obtained from the intensity au- 
tocorrelation functions could be related, via the static 
structure factor S{q), to the hydrodynamic functions pre- 
dicted by a many-body theory derived by Beenakker and 
Mazur (BM) d i]. The study in Ref. [3 showed and 
excellent match between TCDLS measurements and the 
BM predictions in fluid suspensions of relatively low ($ < 
0.35) volume fractions. On the contrary, at higher con- 
centrations ($ > 0.4) the TCDLS data deviated signifi- 
cantly from the BM theory. 
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In the present study we further investigate these phe- 
nomena using a complementary experimental technique, 
that is X-ray Photon Correlation Spectroscopy (XPCS) 
(for some recent reviews, see pi]| - fl2| and references 
therein). XPCS is the equivalent of DLS in the X-ray 
domain. It is not affected by problems related to multi- 
ple scattering and access to larger momentum transfers q 
is possible thanks to the shorter wavelength. An impor- 
tant consequence of the access to higher q values is that 
the static structure factor can be experimentally deter- 
mined and modeled, e.g. by the Percus-Yevick formal- 
ism. XPCS can only be performed at the latest gener- 
ation synchrotron radiation sources where the coherent 
flux is large enough to perform scattering experiments 
and obtain a good signal-to-noise ratio of the correlation 
functions. Damage to the sample induced by the X-ray 
beam can be a nuisance, especially when studying soft- 
matter or biological systems, and during an experiment 
one must carefully monitor the state of the sample. 

Here we report the results of XPCS experiments on 
fluid suspensions of sterically-stabilized spherical parti- 
cles with volume fractions up to $ « 0.49. Previous 
XPCS experiments on hydrodynamic effects in colloidal 
suspensions have focused on charge-stabilized particles 

II. EXPERIMENTAL DETAILS 

The XPCS experiments were performed using par- 
tially coherent X-rays at the IDIOA beamline (Troika) 
of the European Synchrotron Radiation Facility (ESRF) 
in Grenoble, France. A single bounce Si(lll) crystal 
monochromator was used to select 8 keV X-rays, hav- 
ing a relative bandwidth of AA/A « 10"'*. Higher or- 
der light was suppressed by a Si mirror placed in the 
monochromatic beam. A transversely partially coherent 
beam was defined by using a set of high heat-load sec- 
ondary slits placed at 33 m from the undulator source, 
a beryllium compound refractive lens (CRL) unit placed 
at 34 m from the source thereby focusing the beam near 
the sample location, at 46 m, and by a set of high preci- 
sion pinhole slits with highly polished cylindrical edges, 
placed just upstream of the sample, at 45.5 m (see fig- 
ure [Ij. The final beam size selected by the beam-defining 
pinhole slits was 10 x 10 /im^. The parasitic scattering 
from the slits was suppressed by a guard slit placed a few 
cm upstream of the sample. Under these conditions, the 
partial coherent flux on the sample was ~ 10*° ph/s. 

The static scattering from the colloidal suspensions 
was recorded by a charge-coupled device (CCD) with 
22 /xm pixels size located 2.2 m downstream of the sam- 
ple. The dynamic information was obtained with a scin- 
tillation detector (Cyberstar) connected to a multiple-tau 
FLEX01-08D hardware correlator from Correlator. com. 
The detection area was limited to a size corresponding 
to a few speckles by precision slits placed in front of 
the point detector. Typical detector slit settings dur- 




FIG. 1: (Color online) Sketch of the experimental setup for 
XPCS. 

ing these experiments were between 50 x 50 /.im^ and 
100 X 100 /xm^. 

The colloidal suspension was prepared by A. Schoficld 
of the University of Edinburgh and consisted in 
poly(methyl methacrylate) (PMMA) spherical particles 
coated with a thin layer of poly-12-hydroxy steric acid, 
suspended in decalin. A net electric charge on the col- 
loidal particles can typically be produced in low polarity 
solvents b y th e addition of surfactants or charge-control 
additives [l3- However, as prepared, the particles sus- 
pended in decalin are expected to show no detectable 
traces of residual charges (see e.g. [l^) and interact as al- 
most perfect hard spheres. This assumption is supported 
by the XPCS measurements presented below which are 
in good agreement with theories assuming a hard sphere 
interaction potential. The batch solution had a measured 
volume fraction of $ = 0.327. Higher concentration sus- 
pensions were obtained using a centrifugation process, 
while lower concentration were prepared by adding more 
decalin solvent. 

The sample cell consisted in a 1.5 mm diameter Kapton 
tube with wall thickness of « 80 yum. A syringe pump 
purchased from Harvard Apparatus connected through 
Teflon tubing and leak tight fittings purchased from Up- 
church Scientific was used to fill the Kapton tube. This 
experimental setup has the advantage of allowing mea- 
surements under continuous flow as explored previously 
[20| . However, all the measurements reported here were 
performed on stationary, non flowing samples. Our re- 
sults on the dynamics under continuous flow will be de- 
scribed in a subsequent publication. The flow option 
was only used to periodically renew the sample by bring- 
ing fresh particles into the beam. Through many re- 
peated measurements it was found that flowing the sam- 
ples introduces new time scales in the measured corre- 
lation functions [2l|, [2^ that, in particular, can affect 
the long-time decays, even for a long time after the flow 
was stopped. As a consequence, a stabilization time was 
always allowed after renewing the samples or filling the 
flow cell. 

In the high-dilution limit, hydrodynamic or direct in- 
teractions between the particles are negligible and the 
colloids undergo Brownian motion with a diffusion coef- 
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FIG. 2: Measurement of Do — l/rq^ by Dynamic Light Scat- 
tering. The squares sliow the (g- independent) data points 
obtained with the sample (dispersed in decahn) which was 
used in the XPCS experiments. For comparison, data taken 
from an identical sample suspended in cis-decalin (see text) 
is also shown. 



ficicnt Do described by the Stokes-Einstein diffusion re- 
lationship 



Do 



(1) 



6Trr/aH 

where rj is the viscosity and an the hydrodynamic radius. 
Do was measured by DLS on samples with $ < 1%. It 
should be mentioned that with this sample concentration, 
it is impossible to perform XPCS measurements due to 
the weak scattering of X-rays. The DLS measurements 
were performed at several scattering angles and wave- 
lengths of A=532 nm and 633 nm. The time constants 
T were obtained by fitting the correlation functions with 
simple exponentials, and the momentum transfer q was 
calculated from the scattering angle 29 using 

Aim 20 , , 

sm— . (2) 



Here n=1.48 is the index of refraction of the solvent - 
decalin, a mixture of 50/50 cis- and trans-decalin as de- 
termined using the viscosity measurements described be- 
low [2^ [2^. According to Pick's law, the mean square 
displacement of Brownian particles from their position 
at t=0 is < Ax^ >= ^Dot with Do being the diffusion 
coefScicnt. In reciprocal space, this corresponds to a q^^ 
dependence of the correlation time t, and a diffusion co- 
efficient Do = l/{Tq^). As seen if figure [2] this quantity 
is independent of the scattering angle, leading to a dif- 
fusion coefficient for the PMMA particles in decalin of 
£10=9x10^ AVs. 

As a cross-check, Do was measured also for a sec- 
ond sample, consisting of the same PMMA particles sus- 
pended in pure cis-decalin, which has a nominal viscos- 
ity of 3.06 cP (at 24 °C, jgl). The results of the DLS 



TABLE I: Measured Do (DLS) and viscosity 77 (U tube vis- 
cometer) for low concentration suspensions of PMMA par- 
ticles in decalin and cis-decalin (at 24° C), and the hydrody- 
namic radius an calculated from the Stokes-Einstein relation- 
ship 



solvent 


n (cP) 






decalin 


2.6 ±0.1 


(9.0 ±0.2) X 10^ 


931 ± 59 


cis-decalin 


3.2 ±0.1 


(7.35 ± 0.07) X 10^ 


925 ± 38 



measurements, are summarized in Table HI The viscosity 
of the dilute suspensions was directly measured for both 
samples using a U-tube viscometer, hence allowing to de- 
termine au from equation[T] The hydrodynamic radii ajj 
measured in the two different solvents are equal within 
the error bars. 



III. RESULTS 
A. Static properties 

SAXS measurements, performed at the IDIOA beam- 
lines (ESRF), were corrected for background scattering 
contributions from the solvent, the Kapton tubes of the 
sample environment, etc. The static data was fitted us- 
ing the Percus-Yevick (PY) closure [2^ for the structure 
factor S'(g), and a form factor P{q) for spherical particles 
obtained from the fits on the lower concentration suspen- 
sions. The expressions for S{q) and P{q) are calculated 
assuming a particle size polydispersity described by the 
Schultz distribution function |26|. In addition, the re- 
sulting form for the scattered intensity I{q) cx P{q)S{q) 
was convoluted with a Gaussian function describing the 
instrumental resolution. An example of the resulting fits 
for I{q) is shown in figure [3] for the $ = 18.5% sample. 
The same procedure was applied for all samples. The fit- 
ting parameters are the particle radius a=890±12 A with 
a size distribution standard deviation (Ta=89±ll A, the 
volume fraction $ of each individual suspension, and 
an overall multiplicative factor measuring the scattering 
cross section of each sample which is not discussed here. 
On a subset of the samples, additional SAXS measure- 
ments were performed on the SAXS beamline ID02 at 
ESRF. The data measured on the two different instru- 
ments (IDIO and ID02) are in excellent agreement. 

The PY SAXS analysis procedure described above pro- 
vides a fitted form factor for the PMMA particles and a 
fitted structure factor for each of the suspensions. An 
"experimental structure" factor is not directly accessible 
from the data but, assuming that the decoupling approxi- 
mation works well for the relatively monodiperse suspen- 
sion of spherical particles studied here [13] , a good esti- 
mate for it can be obtained by dividing the experimental 
scattered intensity I{q) with the fitted form factor. 



S{q)^I{q)/Piq). 



(3) 
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FIG. 3: (Color online) Example fit for the static scattering 
from the $ =18.5 % sample. The SAXS signal (plotted here 
as I{q)q^ to emphasize the agreement with Porod's law) was 
obtained by circular averaging the 2D CCD images (inset). 
The continuous line shows a fit with a model assuming a poly- 
disperse suspension of uniform spheres interacting via a hard 
sphere repulsive potential, as described in the text. 
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The structure factors S{q) for several samples with dif- 
ferent concentrations can be seen in figure |H where the 
fitted S{q) (continuous fines) are shown together with the 
ones calculated form the experimental data points using 
equation|31 The fitted volume fraction is indicated on the 
graph for each of the individual samples. The agreement 
is very good for all concentrations and over the whole 
q-range. We attribute the small discrepancies that ap- 
pear in some of the fits to experimental artifacts such as 
parasitic scattering. 



qa 

FIG. 4: (Color online) Structure factor S{q) for various sam- 
ples. Data points are calculated, as described in the text, 
from the measured scattered intensity and the polydisperse 
particle form factor resulted from the SAXS fits. Error bars 
are estimated to be smaller than the symbols. The solid lines 
show the fitted structure factor with the PY model for a poly- 
disperse suspension. 



B. Dynamic behavior 

The intensity fluctuation autocorrelation functions 

where measured for wave vectors q around the main peak 
in the structure factor S{q) in a range of 1.5< qa <6. 

Assuming a Gaussian distribution of the temporal fluc- 
tuations at a flxcd q, the normalized dynamic structure 
factor, or intermediate scattering function (ISF), 

S{q,t) 



Siq,0) 



(5) 



is related to the intensity autocorrelation functions via 
the Siegert relationship, 



g^^\q,t)^l + l3 g^'>{q,t) 



(6) 
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FIG. 5: (Color online) Intermediate scattering function 
g^'^\q,t), measured at $ ~ 12.5 % and at $ ~ 48.5 % for 
several q values around the structure peak. Solid lines show 
fits with equation [7] for the low concentration sample and 
fits for the initial decay of the correlation functions measured 
from the high concentration suspension using equation 1101 as 
described in the text. 
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Here, /3 is the optical contrast, which depends on the 
transverse coherence lengths of the X-ray source and the 
sample, and on geometrical parameters such as the pin- 
hole size and the detector slit opening. In the experi- 
ments described here, the contrast /3 was around 2-5%. 

In low concentration samples, the ISF measured at a 
wave vector g is a simple exponential decay, with the 
relaxation rate depending on the single particle diffusion 
coefficient -Dq and q 



(7) 



For higher concentration suspensions, both direct in- 
teractions (DI) acting via the hard-sphere interparticle 
potential and hydrodynamic interactions (HI) between 
colloids mediated by the solvent start playing an increas- 
ingly important role, slowing down the diffusive dynam- 
ics of the particles. While the quasi-instantaneous HI 
between particles are related to the structure factor S{q) 
they do not, themselves, determine the equilibrium static 
structure, nor do they shift the glass transition concen- 
tration (see e.g. [l^). However, the effects of the HI are 
very important, even at relatively low volume fraction as 
they slow down considerably the short-time relaxations. 
Our XPCS results on this complex many-body process 
are shown in the following sections. 



1. Short-time dynamics and Hydrodynamic Interactions 

An important time scale arising in dense colloidal sus- 
pensions is the short time limit ts. This is usually asso- 
ciated with the random motion of individual particles in 
cages formed by neighboring particles. For t > ts the dif- 
fusive motion is slowed down by both HI and DI between 
the particles. At t < ts, the motion is still diffusive but 
faster, as it is slowed down only by HI and not yet by 
DI. This effect can be clearly seen in figure [H The ISFs 
arc plotted here for three different values of q and two 
different samples - a low-concentration one, $=12.5 %, 
and a high concentration one - $=48.5 %. 

While the low concentration sample shows single ex- 
ponential decays, in the high concentration samples 
two distinct relaxation rates are observed. On the 
semi-logarithmic scale used in figure [Sj the exponential 
ISFs g^^''((7,t) appear as straight lines, with the relax- 
ation rates measured by their slopes. A short-time of 
Ts « 2 ms can be estimated from the change in slope ob- 
vious with the high concentration sample (but absent in 
the low concentration one). A physical interpretation of 
this time scale emerging in high density suspensions, was 
given by Segrc, Behrend, and Pusey Q. There, rg is the 
time required for a particle to diffuse away from the posi- 
tion at t~0 to a distance equal to its radius. With our ex- 
perimental parameters, this leads to ts = o^/Dq ~ 9 ms. 
While this is within the same order of magnitude with the 
value observed experimentally, the agreement is clearly 
not very good. An obvious problem associated with this 
estimate for ts is associated with the fact that the free 



diffusion coefficient Dq was used while the diffusion is 
known to be slowed down by HI. Replacing Dq with a 
slower diffusion coefficient (e.g. Ds{q) - see discussion 
below) makes the disagreement between the experimen- 
tally observed ts and the calculated one even stronger. 
This is, however, not surprising, because the length scale 
over which the particles can move before being affected 
by DI should be dependent on the volume fraction, and 
is not necessarily equal to the particle radius. A differ- 
ent estimate of the short time ts can be achieved using 
an expression for the frequency of collisions derived from 
a theory by M. Smoluchowsky which is usually used to 
describe coagulation kinetics (although here we assume 
that coagulation is prevented by the steric stabilization 
of the particles). For a suspension of N colloidal particles 
in a solvent of dynamics viscosity rj and a total volume 
V, the rate of collisions is given by (see for e.g. [2^, |30| 
and references therein). 



V 



V 



(8) 



With a particle radius a, and volume fraction $, this 
results in a time between collisions, identified here with 

TS of 



1 



TS 



(9) 



For $=48.5 %, and using the values of gh and rj from 
Table U equation |9] leads to a short time of Ts ~ 0.5 ms. 
Since even on these short time scales the dynamics is con- 
siderably slowed down by HI the actual time is expected 
to be slower. Considering a factor of «10 for this slow- 
ing down, in agreement with data shown further below 
for the 48.5 % sample, a time scale of rg ss 5 ms can be 
estimated, which is in better agreement with the data in 
Figure [S] 

In the short-time limit, the ISFs can be described in 
terms of a g-dependent diffusion coefficient 



,(1) 



{q,t) = exp [^Ds{q)q^t] 



with the (short-time) diffusion coefficient given by. 



Dsiq) = Do 



H{q) 
S{q)- 



(10) 



(11) 



Here, the hydrodynamic function H{q) describes the HI. 
In the high dilution limit, S{q) = 1 and H{q) = 1, leading 
to Ds{q) = Dq. a non- unitary hydrodynamic function, 
H{q) 1 is a, hallmark of HI. The short-time diffusion 
coefficients were determined by a first cumulant analysis 
[3l| in which the initial decay of the correlation function 
was fitted to an exponential form. The fits were per- 
formed in two stages. First, the correlation functions are 
fitted over the entire time range with a stretched expo- 
nential form. 



(q,t) =/3exph2(rt)^]+(7o 



(12) 
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FIG. 6: (Color online) Hydrodynamic functions H{q) vs qa. 
The data points are extracted from the fitted short-time dif- 
fusion coefficient Ds{q) and the static structure factors S{q), 
using equation [11] The volume fractions for the different sus- 
pensions are indicated on the graph. The solid lines are the- 
oretical predictions of the 5-"f expansion . 

to obtain accurate values for the experimental contrast (3 
and baseline Qoc (with ~ 1 for all the correlation func- 
tions). Subsequently, these parameters are being fixed in 
a second fit, performed only for the initial decays t < tg, 
with a simple exponential form, 

g(^H.q,t)^Pexp[~2Ds{q)qh] + g^, (13) 

where the g-dependent short-time diffusion coefficient 
Ds{q) is the only free parameter. In these ergodic sys- 
tems where the intensity fiuctuations are well described 
by Gaussian statistics, equations [13] and [10] are equiva- 
lent. 

The (5-7 expansion proposed by the BM theory pro- 
vides one of the most successful tools to date describ- 
ing HI in dense but fluid colloidal suspensions. The 
only input parameter required is the static structure fac- 
tor S{q). The BM theory predictions were verified by 
XPCS in several different charge-stabilized suspensions 
with screened electrostatic interactions [l^ [13]. In the 
case of sterically-stabilizcd suspensions, interacting via 
a hard-sphere potential, the predictions of BM theory 
where verified by the two-color DLS experiments de- 
scribed in One problem associated with using light 
scattering is that in general it is impossible to reach high 
enough values of the scattering vector q to obtain accu- 
rate structural information, hence calculated values for 
S{q) (using the PY closure) are usually taken as input. 
In X-ray scattering experiments, access to S{q) down to 



a fraction of the colloidal length scale is straightforward 
so here the hydrodynamic functions H{q) were calcu- 
lated from equation [TT] using measured static and dy- 
namic data. One of the problems encountered in XPCS 
experiments is that accurate dynamic data can only be 
obtained in relatively small g-range {e.g. compared to 
DLS) because the signal-to-noise is often limited by the 
decreasing scattering cross section at high q (or, cquiva- 
lently, by a limited intensity of the coherent X-ray beam) . 
The results for H{q) can be seen in figure [6] for a wide 
range of volume concentrations and a range of q cover- 
ing 2.14x10-31-1 < q <6.44xlO-3A-i (or 2< qa <6). 
Experimental data measured at different concentrations 
are shown by the different symbols specified in the leg- 
end. Continuous lines show predictions of the 5-^ expan- 
sion with no adjustable parameters other than the cor- 
responding volume fractions resulted from the S{q) fits. 
The agreement between theory and experiment is excel- 
lent for solutions with $ < 0.4 which is well in line with 
the earlier observations in Ref. The 5-^ expansion 
employed in the calculation is expected to break down for 
suspensions at high volume fractions [9[. This is indeed 
confirmed by our data, which shows a clear overestimate 
of H{q) at $=48.5 %. 

Figure[7]shows the value of the Hydrodynamic function 
near the structure factor maximum H{qm) as a function 
of the volume fraction. Within the (5-7 expansion, the $- 
dependence of the peak value of H{q„i) for hard spheres 
with PY input for S{q) and $ < 0.45 is well parametrized 
by a quadratic form 

iJ(g„0 = 1 - 2.03$ + 1.74$2, (14) 

represented by the continuous solid line in figure [7] The 
results are also compared with predictions of numerical 
calculations using the Accelerated Stokesian Dynamics 
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FIG. 8: Normalized siiort- and long-time diffusion coeffi- 
cients Ds{qm) / Do, Di{qm)/Do measured at the peak of 5(g) 
{q = q-m) versus the separation parameter |$ — $g|. The con- 
tinuous and dashed straight lines are mode-coupling theory 
predictions as described in the text. 

(ASD) algorithm by Banchio et al. [12, ^ for the HI of 
hard spheres. Their numerical results predict a linear de- 
pendence of H{qm) as a function of $ well approximated 
by 

H{q„,) = 1 - 1.35$, (15) 

which is represented by the dashed straight line in fig- 
ure [71 As it can be seen, the ASD numerical simulations 
provide good predictions over the entire range of concen- 
trations and really excellent predictions above $ >0.4 
where the S-'-f theory fails. This is in agreement with the 
results reported in [l7| for charged-stabilized particles 
with screened electrostatic interactions. 

The dynamic of fluctuations in suspensions with in- 
creasing volume fraction, is also expected the be well 
described by the Mode Coupling Theory (MCT) of the 
glass transition in hard-sphere colloids [sj, [s^. In its 
lowest order in the separation parameter \^g — <I>|, the 
MCT predicts the existence of two divergent time scales 
describing the dynamics in different windows of time - 
the /3-relaxation at fast time scales and the a-relaxation 
at slower time scales. The short-time diffusion coefficient 
Ds{qm) corresponding to the fast /3-relaxation measured 
near q = Qm follows a power scaling law, 

Ds{qm)'x\<S>,-<S>\''\ (16) 

This form is valid both in the liquid, $ < $0 (the sys- 
tems probed here) or glass, $ > <I>g states [3J]. The 
fundamentally different a-process, which is completely 
frozen in the glassy state, restores ergodicity in the high 
concentration liquid phase with a characteristic long-time 
diffusion coefficient Di{q) following 

A(g™)oc|ci>g-<i>|'■^^ (17) 



The slow relaxation process in the high density fluids 
described here is analyzed in more detail in the follow- 
ing section. Here we investigate the scaling forms pro- 
posed by the MCT - equations [HI [T7| - for the diverging 
short-time and long-time relaxations. The results can 
be seen in figure [8] For the low density suspensions, 
the diffusion coefficients obtained from single exponen- 
tial fits arc shown by the open black squares. The short- 
time diffusion coefficients obtained from the high density 
suspensions using the first-cumulant analysis described 
above are shown by the solid black squares. We also 
show additional data (solid black triangles) for the short 
time diffusion coefficient previously obtained by Zontone, 
Moussaid et al. [1^1 in a very similar system - PMMA 
hard sphere particles suspended in czs-decalin. All the 
data follows well the scaling law predicted by the MCT 
and represented by the black continuous line in figure [8] 
with $g = 0.585. 

To date, there are fewer points available for the long- 
time diffusion coefficients. These are shown in figure[5]by 
the grey solid circles, which are also in good agreement 
with the MCT scaling form (grey dashed power law). 
The slow diffusion coefficients were obtained from the fits 
shown in figure [9] and described in the following section. 

The experiments presented in this section provide fur- 
ther evidence for two important results on the dynamics 
in high density colloidal suspensions interacting with a 
hard-sphere potential: 

At moderate concentrations $ <C $c, (i.e. $ < 40 %) 
the short time diffusion coefficients and HI are well de- 
scribed quantitatively by the BM theory. At higher con- 
centrations, only the ASD numerical results by Banchio 
et al. provide accurate results. 

ii) The MCT of the colloidal glass transition provides 
correct accurate quantitative predictions for both the 
short-time and the long-time diffusion coefficients in col- 
loidal suspensions of higher concentrations ($ < <I>g). 



2. Long-time and short-time behavior 

In this section we take advantage of the fact that both 
the short-time and long-time diffusion coefficients are 
readily available from the XPCS data, to test an approx- 
imate scaling law, first proposed by Segre and Pusey [STf . 
By using the aforementioned TCDLS technique, they ev- 
idenced a proportionality between the short- and long- 
time diffusion coefficients measured in concentrated sus- 
pensions with $ «0.46 and higher, over a broad range of 
q values (excluding the smallest qs). This proportionality 
results in a collapse on a single master curve of the entire 
intermediate scattering functions measured at different q 
values when scaled by their short-time decays. This find- 
ing suggests that the structural relaxations of particles 
or "cages of particles" are both related to self-diffusion, 
which contradicts a MCT picture where the a- and /3- 
rclaxations have differe nt p hysical origins. However, sub- 
sequent MCT results [2^ provided a semi-quantitative 
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FIG. 9: (Color online) Intermediate scattering function 
g^'^\q,t) measured at $ ~ 41.3 % (a) and 48.5 % (b) for sev- 
eral q values around the structure peak and fits with double 
exponential decays (continuous lines). 



argument for the approximate collapse of the ISFs on a 
single master curve at high-enough values of q and low- 
enough values of time, where the non-diffusive character 
of the a-relaxations is not very pronounced and the two 
modes hence follow similar decays. The MCT results 
also point to the fact that this scaling form is not ex- 
pected to be valid in charge-stabilized suspensions due 
to the long-range repulsive nature of the potential. This 
was indeed observed experimentally in previous XPCS 
experiments by Lurio, Lumma, Mochrie et al. [isl . . 
Interestingly, a recent study by Martinez et al. [38| using 
XPCS and DLS did not manage to unambigously detect 
the long-time dynamics away from the peak of S{q) in 
high concentration suspensions of PMMA particles, and 
hence could not validate the Segrc-Pusey scaling. How- 
ever, their study, together with an earlier experiment by 



FIG. 10: (Color online) Intermediate scattering functions for 
all values of q measured (same data as in figure [9]) scaled by 
the corresponding short-time relaxation rates Fj, = Dsq^ at 
41.3 % (a) and 48.5 % (b). The insets show the ratio be- 
tween the fitted short-time and long-time diff'usion constants 
versus qa. For both samples this ratio is, within our experi- 
mental accuracy, a g-independent constant. 



Riese et al. [16[ provide a clear experimental proof of the 
equivalence between XPCS and DLS measurements. 

The results presented here are obtained on a system 
very similar to that used b y S egrc and Pusey, or more 
recently by Martinez et al. [38|, namely stcrically stabi- 
lized PMMA uniform spheres (albeit about 1/2 the size), 
and by use of XPCS. As pointed out before, with the 
smaller particles, the higher concentration suspensions 
show clearly two different relaxation time scales within 
the time and length scales accessible in the XPCS exper- 
iments. In order to extract them from the experimental 
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data, the ISFs are fitted with double exponential deeays 
g<''\q,t) ^ Acxpi-T^t) + {I - A)exp{-T2t). (18) 

The fits for the two high concentration suspensions 
studied here with $=41.3 % and $=48.5 % can be seen 
in figure m From this procedure, both Ds{q) and Di{q), 
the short-time and long-time diffusion coefficients, asso- 
ciated with the MCT /3- and a-rclaxations, respectively, 
are readily available. 

In order to test the Segre-Pusey scaling relationship 
between and Di, In [g^^Hq,t)/{Dsq^)'\ are plotted in 
figure lTOl as a function of time like in the original reference 
[37| . These are the exactly the same correlation functions 
as shown in figure [9l except that the above scaling has 
been performed and a log-lin scale is used. As it can be 
seen, the correlation functions measured over more than 
two decades in time collapse on a single master curve, in 
agreement with the scaling proposed by Segre and Pusey. 
The insets in figure [TO] show the ratios between the short- 
and the long-time diffusion coefficients {Dg/Di) at sev- 
eral values of q. As expected according to the putative 
scaling behavior this ratio is a (/-independent constant 
within the experimental error bars. 

We conclude that the Segre-Pusey scaling behavior is 
validated within the combination of length- and time 
scales accessed here, for dense colloidal suspensions with 
$ «40-50 % interacting via a hard-sphere potential. 



IV. CONCLUSIONS 

In summary, we have used a combination of XPCS, 
SAXS and DLS data to measure the (/-dependent short 
time diffusion coefficients Ds{q) and the hydrodynam- 
ics interaction functions H{q) in dense colloidal suspen- 
sions with a hard-sphere interaction potential. Our re- 
sults show good agreement with the BM analytical theory 
at moderate volume fractions ($ <0.4). The XPCS data 
shows good agreement with the ASD numerical scheme 
proposed by Banchio et al. over the entire concentration 
range, and in particular above $ si40 % where the BM 
theory is less accurate. 

The short time diffusion coefficients measured at differ- 
ent volume fractions scale with the separation parameter 
$ — $£, as predicted by the Mode Coupling Theory. 

Finally, the XPCS results for the short-time and long- 
time dynamics in high density suspensions are in good 
agreement with the scaling relationship proposed by 
Segre and Pusey. 

We wish to acknowledge the ESRF ID 10 and ID02 
beamlines for providing the beamtime, T. Narayanan for 
help with the ID02 SAXS experiments and A. Schofield 
for the sample preparation. The work at Brookhaven 
National Laboratory was performed under Contract No. 
DE-AC02-98CH10886 with the U.S. Department of En- 
ergy. 
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